Models of Superconducting Cu:Bi 2 Se 3 : 
single versus two-band description 
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Starting from a model Hamiltonian for the normal state of the topological insulator Bi2Se3, we 
construct a pseudospin basis for the single-particle wavef unctions. Considering weak supercon- 
ducting pairing near the Fermi surface, we express the recently proposed superconducting order 
parameters for Cu doped Bi2Se3 in this basis. For the odd parity states, the d(fc)-vectors specify- 
ing the order parameter can have unusual momentum k dependence for certain parameter regimes. 
Some peculiar results in the literature for surface states are discussed in light of the forms of these 
d(fc)'s. Properties of the even parity states are also illuminated using this pseudospin basis. Results 
from this single-band description are compared with those from the full two-band model. 

PACS numbers: 74.20.-z, 74.20.Rp, 73.20.At 
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I. INTRODUCTION 

The recent prediction 1-4 of the existence of three- 
dimensional topological insulators (TI) and their 
experimental confirmation in Bi2Se3 and related 
compounds^ - — have generated a lot of excitement, as 
these TI form a new class of material which is distinct 
from ordinary band insulators, metals etc, possessing 
peculiar properties such as topologically protected sur- 
face states and usual electrodynamic s 2 ' 9 ' 10 . These prop- 
erties arises from spin-orbit coupling, leading to band- 
inversions at some regions of the Brillouin zone. In- 
terestingly, Bi2Se3, when doped with copper, is found 
to become superconducting— ^—i 1 ^ Some unusual prop- 
erties, such as existence of zero bias conductance peak in 
tunneling experiments 1 ^ 1 ^ and absence of Pauli limiting 
in upper critical field 1 ^, seem to suggest unconventional 
character of the Cooper pairing, though the situation 
is not without controversy^ 9 -. There is a lot of atten- 
tion to the theoretical aspects of superconductivity in 
this compound, in particular possible odd parity pairing 
states. Early on, Fu and Berg2°-, starting from an effec- 
tive Hamiltonian for the normal state of Bi2Se3 near the 
zero momentum T point with two orbitals per unit cell, 
considered various models of superconducting states with 
momentum-independent pairs when expressed in terms 
of this basis. Properties of these models have subse- 
quently analyzed by many, mostly focusing on the surface 
statesi&Z 1 " 2 ^. 

On the other hand, superconductivity in systems with 
strong spin-orbit coupling have been much studied in the 
past, in particular in the context of heavy fermiona 2 ^— . 
There, the usual language used is that the normal quasi- 
particles are described by a pseudospin basis, obeying 
certain symmetry properties, and then the supercon- 
ducting state order parameter and pair wavefunctions 
are expressed in terms of this basis. The superconduct- 
ing states then can be classified by crystal symmetries 
into group representations, and the pairing states are ex- 
pressed in terms the sets of basis functions appropriate 



to the relevant group representations. In this formula- 
tion, for crystals with inversion symmetry, even parity 
superconductors always have pairing wavefunctions even 
in momentum k and are singlets in pseudospins, whereas 
odd parity superconductors have pairing wavefunctions 
that matrices in pseudopsin space with each opponent 
odd in k. This matrix structure is usually expressed in 
terms of a "d- vector" which specifies the corresponding 
pseudospin structure via familiar Pauli 2x2 matrices. 
Properties of the superconductors can then be directly 
obtained by examining these order parameters, includ- 
ing the possibility of surface states. 



Many questions then arise. If the chemical potential 
fi, of the doped Bi2Se3 is sufficiently large compared 
with the pairing potential A, (which is likely to be the 
case since \i ~ OAeV according toi 2 -, whereas the tran- 
sition temperature is < 4K. The measured gap in tun- 
neling is indeed ~ meV, though the precise value is 
controversia l 16 ' 19 ), pairing should effectively only take 
place within one normal state band. How well then 
can we understand the superconducting properties of 
Cu:Bi2Bc3 within a single band picture? What is the 
pairing order parameter in the pseudospin basis? Can 
we understand the Andreev bound states in this way? In 
particular, what is the origin of the very peculiar disper- 
sions found for the odd parity states found i n 16 ' 21 " — ? 



In this paper, we report such an attempt. In Sec |TT1 
we review the model of2P.. Pseudospin wavefunctions are 
constructed in Sec Mil They would be applied to the su- 
perconducting phases, first for the bulk in lPV A| then the 
surface states in HVBl We shall show that many of the 
results in the literature can be understood in this way. 
We conclude in [V] The Appendix gives further discus- 
sions on the surface states and topology of some of the 
superconducting phases. 
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II. MODEL 

In this section, we review the model o&^. We begin 
with the normal state. The effective Hamiltonian which 
captures the physics near k = is given by 

H N (k) = ma x + v z k z a y + va z (k x s v - k y s x ) (1) 

Here fc, k x ,k y ,k z represent the wavevector and its com- 
ponents, v z , v are velocities, a z = ±1 represents the 
two (mainly) p z orbitals in the quintuple layer of Bi2Se3, 
and s x y z are the Pauli matrices for the spin. Equa- 
tion (p} is basically a Dirac Hamiltonian. The ener- 
gies for the quasiparticles are — ±e^ where = 

(to 2 + v 2 k 2 + i^fc 2 ) 1 / 2 with /c|| = ^Jk 2 + A: 2 , thus consist 

of a conduction and a valence band. Bi2Se3 possesses 
D^ d (i?3m) symmetry, which includes parity. Eq ([1} is 
indeed invariant under the parity when this operator is 
taken as oJ&: Eq ([T]) is left unmodified if we substitute 
a VtZ — > —<Jy )Z , k — > —As, with all other variables unaltered. 
Equation ([T]) is actually invariant under the higher sym- 
metry group Dooh as it is obviously unchanged under any 
continuous rotation about the z-axis instead of only 2ir/3 
for D^d- One can check easily that ([lj is invariant under 
reflection about any vertical reflection planes. Since it is 
rotationally symmetric about z, it is sufficient to check 
any one single vertical reflection plane. For example, un- 
der reflection in the x-z plane, k y —k y , s x z — > — s x>z , 
(fTJ) remains indeed unchanged.— 

To discuss the surface states associated with the 
TI, boundary conditions for the wavefunction are 
needed. Unfortunately, this point seems to be somewhat 
controversia l 22 ' 32 ! 33 . For definiteness, we follo w 20 ' 22 here. 
The boundary condition for the wavefunction |$ > at the 
z = plane for a crystal occupying z < is taken to be 
a z \^f >= >, and hence it has no projection in the 
(j z = — 1 orbital. Topological surface states in the form 
of a Dirac cone exist when sgn(mv z ) < O^ 2 - For spin along 
zxk, the bound state energy is Ef, = vkn . Hence, the pos- 
itive energy branch has spin along vz x k. To account for 
the situation of Be2Si3^ we need to take v < 0, though 
we shall consider arbitrary relative signs of v, to and v z 
below for comparison purposes. Since the bulk energies 
are given by ±(to 2 + v 2 k 2 + v 2 k 2 y/ 2 , and so within this 
model, the surface states are always separated from the 
continuum for any given fcii . 

Now we consider the superconducting states, first for 
the bulk. For time-reversal and inversion symmetric sys- 
tems, there is a pair of degenerate states at any given 
momentum k, forming a pseudospin 1/2. Cooper pair- 
ing occurs between opposite momenta k and —k, and 
can be classified into even parity, pseudospin singlet and 
odd parity, pseudospin triplet states^ These supercon- 
ducting states can further be classified by their different 
symmetries under the crystal symmetries into different 
representations in group theory^ - — These representa- 
tions depend only on the point group (but not the space 



group). Possible forms of the corresponding momentum 
and pseudospin dependence in each group representation 
expressed in the form of basis functions: the general form 
of the order parameter can be a linear combination of 
the independent basis functions of the same symmetry, 
each term possibly multiplied by a momentum dependent 
function which is invariant under the particular group 
under consideration. They have in particular been listed 
for the cubic Oh, tetragonal D^h, and hexagonal Dq^ 
groups 2 ^— . For D 6 h we have the group representations 
Ai, B\, A2, B2, Ei and E%, with each of the above ei- 
ther even (g) or odd (u) parity. The corresponding table 
for D^d appropriate for Bi2Se3 was not listed in these 
references, but can be trivially obtained from those for 
D§h since D§h would reduce to D^d if we discard rota- 
tions about I of odd multiples of 2ir/6 and three of the 
horizontal rotational axes. In this case, A\ is no longer 
distinguishable from B\, and similarly for A2 and B2, 
and Ei and E% ■ The resulting group representations and 
their basis function arc listed in the first two columns of 
table QJ following Ref 29 . For simplicity, we do not list 
all the possible independent basis functions, but mainly 
those which would appear again in the later part of this 
paper. For the complete basis function set, we refer the 
readers to the literatur o 29 ' 30 . 

In 2 ^, various types of momentum independent (local) 
pairing in the orbital and spin basis of eq ([U were con- 
sidered. Since this formulation involves states of eq ([T]) 
without a priori distinguishing the conduction and va- 
lence band, we shall refer to this as the "full two-band 
description" 34 The symmetry of these states were already 
discussed in2£, and we list each these states with their 
corresponding symmetries in column (iii) of Table |H In 
this table, we have followed the notation of2°- and use 
1,2 to label the two orbitals instead of <j z = 1,-1 of 
eq (fT|). For convenience of comparison with other works 
in the literatur e 16 ' 21 ' 23 , we also list in columns (iv) and 
(v) these order parameters in matrix form, which we 
shall denote A 1 and A 11 and referred to as "Nambu 
I" and "Nambu II". A 1 is the matrix order param- 
eter in the ordinary Nambu notation after generaliza- 
tion to two orbitals, that is, if we use the operators as 
(c a f, c^.),, cL,, c' a 1 ), where c CTS and cj. s 's are annihilation 
and creation operators, and a = ±1 the two orbitals. 
If we use instead (c^, Co.j,, c CT ,, — cL.), as done in^ 2 -, (or 
(<Vf , c CT 4, — j,,c^), as in 2 ^) the order parameter ma- 
trix is A . The two notations are related simply by 
A 1 = A u (is v ) ( A 1 = A n (-iSy) ). The factorizing out 
of (is y ) to the right has the effect of what has been done 
in the 3 He literature 35 , where the order parameter matrix 
is written as (d ■ s)(is y ), so that d transforms as a vector 
under spin-rotations. In this way, it is clear from column 
(v) that the two entries listed under E u are related by a 
7r/2 rotations about the z-axis. In making this table, we 
have made use of the gauge symmetry of superconduc- 
tivity to simplify the matrices (by removing factors like 
±1 or ±i). However, we have kept the correct relative 
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(i) 


(ii) 


(iii) 


(iv) 


(v) 


even parity 




1 


[1 1 1 4> +[2 1 2 4> 


ISy 


1 






|1 1 2 4-> — |x 4. 2 t> 


%(J x Sy 




Mg 


Im k% 








E a 


Re k z k+ 










Im 








odd parity 


Am 


k z z\ k x x -\- ky*y 


|lt2|> +|H2t> 




cr y s z 


A 2u 


k x y kyX 


ltU> - 2 t 2 4-> 


OzSy 




E u 


Re k+z; k z r+; k\k z r- 


1 2 1> -|im>) 


— <?yS z 






Im 


|1 1 2 t> +|H24.> 


Wy 


CySy 



TABLE I: Representations (column (i)), basis functions (column (ii)) and pairing wavefunctions (column (iii)) for the super- 
conducting phases considered in this papei^. The matrix form of the order parameters for (iii) are given in (iv) and (v) for 
the Nambu-I and Nambu-II representations. Here f± = x± iy, k± — k x ± iky. 



have the correct relative transformation properties. 
I 

I 



phase between the two partners within E u , so that they 



Table Q] shows the correspondence between each local 
pair and the possible basis functions. However, it does 
not tell us directly what exactly the momentum depen- 
dences are, since multiplication of any basis function by 
a function invariant under the crystal symmetry is also 
an as good basis function. It also does not tell us, in 
particular for the case of odd-parity pairing, what lin- 
ear combinations between the different inequivalent ba- 
sis functions (e.g. k x x + k y y and k z z for A\ u ) that we 
should take. We shall see later that this information is in 
fact useful for understanding the surface state spectrum. 
To do this, we have to first construct the single parti- 
cle pseudospin wavefunctions obeying the correct crystal 
symmetries, as we shall do in Section Hill We shall then 
take up the task of showing how the pairings listed in 
column (iii) correspond to the basis functions listed in 
column (ii) in Sec IIV Al 

III. PSEUDOSPIN WAVEFUNCTIONS 

Now we construct the pseudospin wavefunctions. For 
each k, we shall denote the two degenerate states by 
\k,a > and \k,f3 >. We shall demand that the corre- 
sponding pairs at — k are related to those at k by the 
relations 

| - k,a> = P\k,a> (2) 
\-k,/3> = T\k,a> (3) 

where P denotes the parity and T the time-reversal. P 
was already taken to be a x , and we shall take T as (— is y ) 
times the complex conjugate. Correspondingly, we have 

\k,P >= PT\k,a > (4) 



and 

\-k,P>=P\k,p> (5) 

Note that T|fc, (3 >= -\-k,a >, asT 2 = -1. \k,a > and 
\k,/3> are also required to satisfy certain rotational sym- 
metry properties, to be treated in details below. These 
requirements basically enable us to roughly think of a 
and (3 as "spin-up" and "spin-down" respectively^ We 
shall first deal with eq © and ([3]), ignoring the rota- 
tional properties for the moment. We shall thus first find 
an intermediate basis \k,a' > and \k,/3' > obeying equa- 
tions @ and ([3]) without worrying about the rotational 
properties. 

For this purpose, let us first introduce spin wavefunc- 
tions which diagonalize the spin part of eq ([T]), using 

\s = zxk>=l=( ) (6) 

\§ = -2 x k >= -L f ie ~^ ) (7 ) 

for spins along ±z x k. Here <j)^ is the azimuthal angle 
of k in the x — y plane. These spin wavefunctions satisfy 
(k x s y — k y s x )\s = ±z x k >= ±fc|||s = ±z x k >. It 
is then straight-forward to diagonalize eq (J]). We shall 
define \k,a' > for k z > (the "northern hemisphere") to 
be the one associated with spin along £ x k: 

\k, a' >^ -Le^ ( \ + Vk * (8) 
1 y/2N \m + iv z k z J \ ie^k J > 

where the first column matrix denotes the part in orbital 
space and the second part denotes the spin space. Here 
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Et is the energy of the particle (which can be ±e^r), and 
N = [2Et(E^ + vk^)] 1 / 2 is a renormalization factor. The 

other state \k, ft' > for k z > 0, as well as states for k z < 0, 
are obtained by the symmetry requirements ©,© and 
dU>22. We thus have, for k z > 0, 



under parity and time-reversal) 

m '* ( K - ( " sh % IAjIcob^W^s) 
x { ' \ |Ag|coB0ge < Ws +o s) sin^g 



(13) 



\k,ft' >e 



ik-r 



V2N 



iv z k z 
+ vk\\ 



and 



— ik-r 



V2N 



TO + ZVzfcj; 

E- k + v fc„ 



le 



-i<t>t 



(9) 



(10) 



m' eff (fc) = 



and 



COSfpr 

|A £ |sin^e l W£ +Q ^ 



|yl £ |sin^e- i ( <#, E +a E) 



m/ ff (fe) = ( _ ieiWs+os) ) (15) 



(14) 



k,ft' >= 



( E k- 


f«fc|| \ 




where 


I TO — 









We have written the last two wavefunctions for wavevec- 
tors in the "southern hemisphere" using the labels — k 
with k z > 0. This is for convenience later since we shall 
always be consider Cooper pairs between k and —k, and 
it is sufficient to write these pairs with k z > for the 
pair due to the fermionic antisymmetry of wavefunctions. 
Note that 



'-k 



We now proceed to find the wavefunctions \k, a > and 
\k,ft > with the desired P, T, and rotational properties. 
One can of course directly study the wavefunction them- 
selves. However, a more convenient way to proceed is to 
evaluate some physical quantity with known transforma- 
tion properties (c./i^l). For this, we consider the spin op- 
erators projected onto our two-dimensional Hilbert space 
\k,a> and \k, ft > for each k point. These operators are 
thus then also 2x2 matrices. This spin operator is related 
to the effective magnetic moment of our quasiparticles, 
if the orbital contributions can be ignored (which can in- 
deed be the case if the relevant orbitals are just p z and 
the mixing to p x ± ip y can be ignored) . We shall there- 



fore denote them as 



Anyway, the operators for this 



effective spin moment m| , j = x, y, z, are simply 

'off /A _ ( < k,a'\sj\k,a' > < k,a'\sj\k, ft' > 
mj [ } ~ \ <k,ft'\s j \k,a' > <k,i3'\ Sj \k,/3' > 

,( 12 ) 

where the Sj inside the matrices are the spin Pauli ma- 
trices as in eq ([1]). The prime ' is to remind us that we 
are using the \k,a' > and \k,ft' > basis at this moment. 
Viewed as operators, we thus have 

m° s {k)= |A J 7><*»7l*j|*.7 / ><*»7 , | 

7,7'— a', /3' 

Straight-forward calculations using eqns flSJ) and © give, 
for k z > 0, (those with k z < can be found later by us- 
ing eq ([H and ([3]) , this guarantees the correct properties 



A n = \A % \e- ia * = —{E % + vk\\)(m - iv z k z ) (16) 

is a factor generated by the overlap of the orbital wave- 
functions in eq ([8]) and (O, and so 



(TO 2 +t-2fc2)l/2 



and 



e~™z = (sgnEj) 



m — iv z k z 
(m 2 + u^fc 2 ) 1 / 2 



(17) 



(18) 



The set of 2 x 2 matrices p\ 2 3 ,— with p\{k) = \k, a' X 

k,ft'\ + \k,f3' >< k,a'\ etc do not yet have the desired 
transformation properties under rotation. We now con- 
struct a new basis \k, a >, \k, ft > so that the correspond- 
ing Pauli matrices p x ,y,z for the pseudospin do transform 
like an axial vector. Since the system has complete ro- 
tational symmetries about z, we must require m c z s oc p z . 
To do this, we simply have to find a basis so that m z is 
diagonalized. This can be done by choosing 

|jfe, a >= ^ a' > -ie^+^&ft' >) (19) 



and 



\k,ft>= 



V2 



(\k,ft' 



a ^\k,a'>^ (20) 



where the phase factor is at this time arbitrary. Note 

that we have demanded that \k, ft > be related to \k, a > 
by eq ((4]). In this new basis, we find 



mf(k) = \A % \p 



(21) 



and 



5 



(|^|cos^ - isin<^)e^E +a E +2f ^ 



.off 







(l^lsin^g + icos0 s )e 



(l^lsin^g - icos^)c 



(23) 



To proceed further, it is simplest to examine the radial and azimuthal components of rh and p, i.e. m 

off i n^nA^meff 

r 



coscf)^m x s + sin0gm® ff and m^ ff = 



-Mn^mj 11 + cos0£m° , and similarly for p r and 



i.e. 



iff — 



and 



le 











Evidently due to the existence of vertical reflection 
planes at arbitrary angles with respect to the x-axis, m* ff 
( m^ ff ) must simply be proportional to p r ( p^) but would 
not involve the other component. One sees that we can 
choose 8j: to satisfy 



2% = 



(24) 



or aj: + = 7r. We shall adopt the first choice. In this 



case we get 



and 



m 



cff 



\Mpt 



.cff 



(25) 



(26) 



For this choice, a pscudospin along the positive azimuthal 
direction would correspond also to an effective magnetic 
moment and hence spin along the same direction. (The 
alternate choice would give m r = — \Aj:\p r and m*^ = 
—ptj, instead.) Back to the Cartesian form, we have 



off 



(fc) = P X - (1 - \A^\)k X {k X p X + kyPy) 



.off 



(k) = fly - (1 - |^| 



kyPy) 



(27) 



(28) 



which explicitly shows that m xy has the same transfor- 
mation properties as p x . y respectively. 

Note that the procedure above also gives us the ef- 
fective <?-factor for the effective moments. For mag- 
netic moment along z and the radial component r, eq 
(|2"Tj) and (|25|) show that they are reduced by the fac- 
tor \A%\ = (m 2 + v 2 z kl) 1/2 /\E^\ < 1 given in eq (fTT]) . 

but there is no reduction for the (f> component. For k 
on the Fermi surface, = p, \A?\ = 1 for k paral- 
lel or antiparallel to z. It decreases for increasing ku, 



and for k in the x-y plane, — > \m/p\, which can be 
substantially less than unity, as in the case relevant to 
the experiments' 12 . This effective moment would be rel- 
evant when considering questions such as Pauli limiting 
of upper critical fielcU^, or spin susceptibilities measured 
by Knight shifts. Returning to the pscudospin basis, eq 
COI and (HOI) become 



l/c, a >= 



and 



-iag/2 



(\k,a' > -ie i( -^ +a ^\k,p' >) (29) 



\k,/3 >= e —^r- (\k,/3' > -ie-^+^l^a' >) (30) 



States at — k, k z > 0, can be obtained by using eq ([2]): 



— k, a > = 



and 



|-fc,/3>= 



e -i« S /2 



V2 



(\ -k,a'> -ie i (+*+ a K)\ - k,f3' >) 

(31) 



"7T 



(\-k,p'> 



-ie 



, a > 



(32) 

With \±k, a' > and \±k,P' > available in eq © © , CEU]) 
this completes our construction of the pseudospin 
basis. We shall express the Cooper pair wavef unctions in 
terms of it in the next section. 

Before we proceed, since we would also be interested in 
surface bound states in the superconducting states in Sec 
IIVB1 we consider reflection of quasiparticles at a surface 
in the normal state before we end this section. We con- 
sider a crystal occupying z < 0, with a surface at z = 0. 
Consider incident wavevector fc = k x x+k y y+k z z, k z > 0. 
Due to our ways of writing wavefunctions for wavevectors 
in the southern hemisphere, it is convenient to write the 
reflected wavevector as — k' where k' — —k x x — k y y + k z z 

and note that 6 p, = (br- 



and use eq (1311) and 
Straight- forward algebra shows that \k,a > is reflected 
only into | — k',a >, and similarly for a — ¥ j3. Indeed, 
the wavefunctions >= \k,a > 
\k,0 



,a > or 



fc',/3 >, with the reflection coefficient 



6 



-(sgngg) , m+iv * k * 



fc; (m 2 +t ,2 fc 2)i/2 satisfy the bound- 
ary condition < a z = — 1 1 \1/ >= 0, since one can easily 
verify that < a z = — l\k, a >— e~ m E < u z = — l\k',a > 
and similarly with a — > j3. Hence the reflection of quasi- 
particles at z = in the normal state does not alter the 
pseudospin species, nor is the phase shift dependent on 
the incident species. Hence, in the model of eq (P), the 
z = surface is pseudospin-inactive, a result that we 
would use in Sec IIVBI 

Our single band description here is unable to capture 
the surface states of a TI. These states are superposition 
of states from both the conduction and valence bands. 
The implication of this for the surface states of the su- 
perconducting phases would be discussed in Sec I1VB1 



IV. SUPERCONDUCTING STATES 
A. Bulk 

It is straight-forward to obtain the order parameters 
for the superconducting states in our pseudospin basis. 
We discuss each of the phases listed in Table Q] in turn. 
We confine ourselves to momentum independent pairing 
within the a, s basis of eq ([I}. Generalization to addi- 
tional momentum dependence is straight-forward. So far 
for Cu:Bi2Se3, superconductivity has been found only for 
H > 0fi£ but we shall also consider general sign of /i in 
the following. 
A lg : 

Both the "intra-orbital opposite spin pairing" |1 j" 
, 1 \> +|2 j, 2 \> and the "inter-orbital singlet pairing" 
|1 1)2 \> — 11 1, 2 |> has Aig symmetry. In general, they 
are expected to be mixed. However, since they have also 
been discussed separately in, e.g.^, we shall also first do 
the same likewise, and consider a general linear combi- 
nation later. For ease of referral, we shall refer these two 
states as A' lg and A lg respectively. 



A' ■ 

The pair wavef unction is |1 1 
,2 l>. The corresponding form 
dospin language is just S'fc :7 , 7 



_i> +|2 t 
in the pseu- 

Here the 



< ky, -fcV|l t, 1 i> + < *7, -*y|2 t. 2 i> 

prime over the sum means that k is restricted to the 
"upper hemisphere" (as those in the other hemisphere 
is already included by antisymmetry), and < 1 j, 1 

-ky' >=< 1 t |&7 >< 1 i | - hi > - < 1 I 
|&7 >< 1 t | - hy' > can be evaluated using 1(29 ]) - ([32 ]) 
and eq (f5])- (fTT]) . We find that this simply reduces to 

J^g \ka, —k/3 > — \k/3, —ka > , with no additional mo- 
mentum dependent factors. If the pairing term in the 
superconducting Hamiltonian is taken as A^(eLcj, + 

cj-scli )+h.c. where c\ s are the creation operators of orbit 
a and spin s and h.c. indicates the Hermitian conjugate, 
then the corresponding term in the pseudopin basis sim- 



ply reads VrA! (ct C* - — etc 1 - ) + h.c This has 

v J IV ]~a -k/3 k/3 -ka' 

just the familiar form for momentum independent con- 
ventional s-wave pairing. The quasiparticle energies in 
the superconducting state, measured with respect to the 
chemical potential /i, are ±Es with E$ in the familiar 
form, for fi > 



El 



(e g - M ) 2 + (Ai) 2 



(33) 



where we have taken the gauge where Ai is real. The 
corresponding formula for /i < is 



E% = (e % + nf + {A' x f 



(34) 



Note that since we started with a normal metal and 
then introduce the superconducting pairing, we necessar- 
ily have /i 2 > to 2 implicitly, and weak-superconducting 
pairing actually requires further that — \m\ 3> lA'J. 
For the full two-band description, the Hamiltonian in the 
Nambu-II notation is just (Hjy— ^)r 2 +A / 1 r a; (see Table[T|, 
where T XtVtZ are the Pauli matrices in particle-hole space. 
Now there are instead two pairs of allowed E$ due to the 
presence of two bands, given by Eg — (eg ± /j,) 2 + (A^) 2 , 
thus including both eq ([33]) and (J34J) irrespective of the 
sign of fx. 

The pair wavefunction is |1 f, 2 \> -|1 |, 2 f>- Fol- 
lowing the same procedure described above gives the re- 
sult J^'g -g^dfca, — fc/3 > —\k(3,—ka >) in the pseudospin 

k 

basis. If the pairing term in the Hamiltonian is written as 
Aj (c{^C2j_ — c\^c*2i)+h.c. with A 1 real, the corresponding 
quasiparticle energies are ±Es with 



(35) 



with the upper and lower signs for fi > and /i < 
respectively. Again the energy gap is isotropic in k space 
and is given here simply by | — A x \ . Note that the system 

becomes gapless if m = even for finite A x . On the other 
hand, in the full two-band description, with Hamiltonian 
Hg = (Hn — (Jt)T z + A 1 <j x t x , we obtain again two pairs 
of energies, given by 



1/2 



Considering the lower energy branch (where eg 



(36) 
±H for 



/x^O) and taking the weak-pairing |/z| — |m| 3> |Ai| ap- 
proximation, we recover eq (133]) . as expected. Generally, 
the system is gapped whenever m ^ and A : ^0. If 
A : = 0, we recover the normal state, and states at mo- 

±(j, have Eg — 0. If m = 0, 
| = [i 2 + (Aj ) 2 can be satisfied, 
and at positions slightly different from the one-band re- 
sult eg = ±/i with correction due to finite A x . 
general A\ g : 



mentum k such that eg 
gaplessness occurs if e| 
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The most general Ai g pairing wavefunction is a linear 
combination of that of A' lg and A 1 . The general pairing 

Hamiltonian is A^c^c^+c^cJjjJ + A-l (cj^c^ — c|^C2*) + 
h.c The corresponding expression in pseudospin lan- 
guage is again just the linear combination of those given 
above for A[ and A l We shall, for simplicity, restrict 



ourselves only to the case where time-reversal symmetry 
is preserved, thus A[ and A x can be chosen to be real si- 



multaneously, though either can be positive or negative. 
The energy spectrum is, in the single band description 



E 2 s = {^ k Tnf 



a; + -a, 



(37) 



The system is gapped if A[ + ^Aj 7^ 0. The correspond- 
ing result in the full two-band description is 



J 



(A' 1 ) 2 + (A' 1 ') 2 ±2 



2 2 
H to 



(A'^) 2 - 2/imAiA; + (e| - m 2 )(^ 2 + K 2 ) 



1/2 



(38) 



r 



The lower energy branch reduces to eq Q37[) in the weak- 
pairing limit. If the full expression (|38[) is used, one can 
check that the system is gapped whenever mA, + A iy u ^ 
0. (Thus recovering the one band result since there /i 
must be finite). If mA 1 + A x fi = 0, gaplessness still 
requires the condition (v ■ k) 2 = (/1 2 — to 2 ) — (A 2 — A^) 
be satisfied. Hence we have the following : 

(1): If fj, ^ 0, we need A^ + 2 A' x ' = and (v ■ k) 2 = 

(/i 2 — to 2 ) ^1 — (A 1 //i) 2 ^j for gaplessness. Note that for 

/i 2 ^to 2 , the latter happens if and only if /x 2 ^A 1 2 . We 



reproduce the weak superconducting pairing results when 
- |m| > Ai and A'i'. 

(2) For /i = 0, the system is gapped whenever to 7^ 
and Ai 7^ 0. (2a) If to = 0, gaplessness occurs only if 
I A'i| > |Ai|. (2b) If A'i = 0, the system reduces to A' lg , 
and the system is fully gapped unless A^ also vanishes. 
We shall use these results when we discuss the surface 
states and topology in Sec IIVBI and the Appendix. 



Aiu 

becomes 



The pairing wavefunction, |1 f, 2 i> +|1 I, 2 f> 
les 



(m 2 + v 2 k 2 y/ 



1/2 X! 



v(k x ik y )^ a _% a > + v ^+ lk v\ kl3~kp > +(sgn^)^(|fca,-fc/3 > +\k(3,~ka >) 



r 



IC-. c - 



If the pairing term in the Hamiltonian is written as 
Ai«( c i-f c 2i+ c ii c 2t)"'~^ -C "> with Aiu real and positive, the 
corresponding form in pseudospin is (— d x +id 

dJcL c f - + ct c f - 

zy ka -k/3 k/3 -ka 

quasiparticles are taken at the Fermi energy ji, 

vk T 



ka — kc 

) + {d x + idy^lpC^p with, when 



\H\ (m 2 + v 2 k 2 ) 1 / 2 



d z = A iu (sgn^) 



v z k z 



(to 2 + v 2 k 2 y/ 2 



(39) 



(40) 



As compared with the Balian and Werthamer (BW) 
stated where d(k)\\k (or d(k) anti-parrallel to k, after a 
gauge transformation), the d{k) here has a very peculiar 
form. The ratio between the in-plane and z component 
is 



Besides the anisotropy factors from the velocities, an ex- 
tra factor m/n arises, which suppresses the d|| component 
relative to d z if |m| < \fx\. This, in retrospect, is actu- 
ally not surprising since the pairing is between opposite 
spins in the original a and s basis. A pure opposite pseu- 
dospin pairing would have d parallel or antiparallel to z. 
The x, y components of d are actually generated by spin- 
orbit coupling. Moreover, we note that the relative signs 
between d XyV and d z depends on the signs of the var- 
ious parameters of the system, in particular sgn(i; 2 TO). 
We shall come back to this when we discuss the surface 
bound states. We note here also that this peculiar rela- 
tive sign and magnitudes between the components of d is 
allowed here due to the inequivalence between z and x-y 
under the relevant D^d symmetry. For a cubic system 
such as YPtBi, 43 k x x + k y y and k z i necessarily comes in 
the combination k x x + k y y + k z z. 



l \\ 



vkn 



(j, v z k z 



(41) 



Despite the peculiar form for d, the energy gap turns 
out to be isotropic in the weak-coupling limit. The square 
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of this gap is given by d • d, which is 



the corresponding d(k) is, in the weak pairing limit, 



A 2 



l> 2 fc 2 



mvku 



+ vlk 2 



which works out to be simply A 2 U (1 — (m//i) 2 ) when 
we restrict ourselves to particles near the Fermi surface, 



where m 2 + v 2 k 2 + v 2 k 2 



For the energy gap, 



the anisotropics due to eq pTj) and the overall factors 
(m 2 +w 2 /c 2 ) 1/2 in eq ((39]) and (gO]) cancel each other. The 
quasiparticle energies are just (eg=F/x) 2 +A 2 M (l — (m//i) 2 ). 
This phase is fully gapped (provided /i 2 > to 2 , which, as 
mentioned, is necessarily the case for a single-band weak- 
pairing superconductivity picture to be meaningful). 

In the full two-band description, the quasiparticle en- 
ergies are already worked out 



4 



TO 



1/2 



(42) 



which reduces to what has been just given in the weak- 
pairing limit. In this two-band description, the sys- 
tem canstill be gapless when A lu ^ provided to = 
±^//i 2 + A 2 U (with gapless point at k = 0). We shall use 
this result later in the Appendix. 
A 2u : 

The pair wavefunction |1 f, 1 l> — |2 t, 2 l> becomes 



v ( k x iky ) 



Er 



ka. —ka > 



v(k x +ik y ) 



\kP,-kp > 



the pairing term is written as A2„(c}^.c^ 



' C 2-t C 2i) 



If 

- /i.e., 



d(k) = A 2 „ 



v(k x y - k y x) 



l< 



(43) 



The magnitude of the gap is just A2 tl |vfc||//^|. This is the 
usual planar phase in the 3 He literature, and is regaining 
attention due to its analogy with topological insulators 
in two-dimensions (e.g . 44 ' 47 ). In three-dimension how- 
ever, this state has point nodes in the gap at the north 
and south poles of the Fermi surface, where k x , k y both 
vanish. 

The expression for the quasiparticle energies in the 
two-band description is 



El 



4- 



1/2 



(44) 

which reduces to the above results in the weak-pairing 
limit. The state is gapped at all vk\\ ^ 0. Gaplessness 



A 2 



i 2 can be 



can occur if the condition v 2 k 2 
satisfied. 
E u : 

This is a two-dimensional representation, as |1 \. 2 ^> 
is the time-reversed of |1 j" 2 / | s >. Generally, the super- 
conducting state can be a superposition of the two. Let 
us first consider the state i\l t, 2 f> ( we have inserted 
an i factor for later convenience.) If the pairing term in 
the Hamiltonian is given by iA u (c[^c^) + h.c, then we 
have 



M) = A ^4 (Sgn ^( (m 2 + ! 2 fc 2 )V 2 



vk + 



(to 2 + v 2 z k 2 z y/ 2 2|^| 



( (ma+ ^2)1/2 " P ) M^)*" 



r 



(45) 



This d(k) is complex (dx d* ^ 0) reflecting the fact that 
the state |1 y 2 t> has broken time-reversal symmetry. 
The first two terms in eq (|45l) are proportional to k z f + 
and k + z listed under E u in Table [J] The last term has a 
more complicated momentum dependence, but since 



f 



1 



{m 2 + v 2 k 2 y/ 2 | m | 



\v 2 B 



for small fcii, it is simply proportional to k z k 2 + r-, the 
third independent basis function listed in Table HI in 
this limit. The spectrum for this state is complicated 



since it is " non- unitary" , that is, the energy of the 
two pseudospin-species at the same k point are typi- 
cally unequal, due to the lack of time-reversal symme- 
try. We shall not investigate this phase in detail, but 
turn to the time-reversal symmetric states within this 
two-dimensional manifold. 

Let us consider then i\\ y, 2 y> |, 2 \>. This 

state is just the linear combination of the one discussed 
above and its time-reversal conjugate. The d(k) vector 
for this state is therefore simply twice the real part of eq 
(1431) . and so 
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;( s S n M) 



(m 2 + 1)2^2)1/2 |^| 



v z k z (cos(20£)i + sin(2<^g)y) 



(46) 



r 



corresponding to the basis functions listed in the first line 
under E u in Table HI Despite its complicated form, the 
square of the gap, obtained from d- d, is given simply by 

2 f 2 fc2 + v 2 z k 2 z 
p 2 

This state has two point nodes, as it is gapless for k 
parallel to y. The result is in accordance with the full 
two-band result, which is 



El = el + f 



■Af±2 



1/2 



(47) 



and is just eq f|44[) with v z k z — > vk y . The gap-squared in 
the weak-coupling limit is 



v 2 kl 



_ A 2 V 2 k 2 x + v 2 z k 2 

,i 2 



for momenta on the Fermi surface. The point node for 
this phase has also been noted in 2 ^. 

d(k) for the state |1 f 2 t> +|1 I 2 I> is evidently 
twice the imaginary part of eq (|45l) . It is just eq (|46|) 
rotated by ir/2 about z. 



B. surface states 

Now we consider the surface states in the supercon- 
ducting phases. We shall focus on the odd parity states 
since they have received more attention in the literature, 
(see however the Appendix) 

For weak superconductors (pairing potential much 
smaller than fermi energy) , surface bound states are most 
conveniently discussed quasiclassically. Bound states can 
be formed at the surface since the quasiparticle with 
incident wavevector ki n sees a different order parame- 
ter from when it is reflected into k out . Since the sur- 
face z = is pseudospin inactive^, the problem maps 
to the evaluation of the quasiparticle bound states at a 
one-dimensional junction, where the order parameter for 
z' < is different from that for z' > 0. Here, the order 
parameter for z' < can be identified with that of k- m , 
and z' > with fc out . In all situations relevant to us, the 
magnitude of the order parameter |A'| for k- ln and fc ou t 
are identical, only the phase £'s are different. The effec- 
tive phase difference for the junction is x — (out — (in- 
The absolute value of the bound state energy is just 
\Eb\ = |A'|| cos(x/2)| if assume that the order parameter 
is constant up to the surface. Since we shall be inter- 
ested also in the sign of Eb, we give a short discussion 
of it££. For normal state with particle-like dispersion, 



where the energy is increasing with the magnitude of the 
wavevector, the bound state energy for positive z' mo- 
mentum is E h = — |A'| cos(x/2) < for < x < 2-7T, and 
E b = |A'| cos(x/2) > for — 2tt < x < 0, and vice versa 
for negative z' momentum. The above signs should be 
reversed for hole-like normal state dispersions^. For our 
normal state, the spectrum is particle-like for fi > 0, and 
hole-like if \i < 0. 

For time-reversal symmetric odd parity superconduc- 
tors, d(k) can be chosen real. It is best to work with the 
quantization axes z' s (no relation to z' above) which is 
perpendicular to both di n = d(k- m ) and d ou t = d(k out ). 
We shall choose z' s to be parallel to e?j n x d ut- Then the 
order parameter with z' s as the quantization axis is given 
by 



'0 d x 







idyi 



hence diagonal in pseudospin space. Consider first the 
"down" component. We can write (d x > s + idy^)^/ out — 

|A'|e Ci »/°ut where Cm/out is J ust tne an S le for ^n/out in 
the x s - y s plane, measured counterclockwise from the 
Hence the phase difference % for the "junction" 
is just x — (out ^ (in = (j the angle between di n and 
cT out , with < ( < 7T, and with E b = -|A'| cos((/2) < 
for particle-like normal state spectrum. For pseudospin 
"up" along z' s , we can rewrite (— d x ^ + id y /)i n / ou t = 



IA' 



and the effective phase difference is x 



(-(out) — (— (in) = — (■ The bound state energy is 
Eb = |A'|cos((/2) > for a particle- like normal state 
spectrum. 

Summarizing, the bound state energy is positive (neg- 
ative) if the pseudospin is parallel (antiparallel) to 
(sgn/i)dj n x G? ou t. We are particularly interested in com- 
paring this sign with the bound state in the normal phase 
of our TI. We recall that, in our model, the energy is pos- 
itive if the spin is parallel vn x k, where n is the surface 
normal (pointing outward from sample). Hence, if we 
focus on the relative sign between the superconducting 
and the normal state, we can state that: 
The bound state dispersion for the normal and super- 
conducting phases has the same sign if (sgn/i)dj n x d out 
is parallel to vn x k. 

We shall call this situation as "regular relative to nor- 
mal" (RN). Conversely, we shall call it "anomalous rela- 
tive to normal" (AN). 

Actually, another meaningful comparison would be to 
the surface bound state for the BW phase where d is 
parallel to k. In that case, then the dispersion has the 
same (opposite) sign as the BW phase if d m x d out is 
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parallel (antiparallel) to h m x /c out . Since k m and fc out 
differ only by the component along the surface normal, 
we see that k in x fc out is simply parallel to nxk in . We shall 
mainly be focusing on the first comparison, though the 
comparison with the BW phase can be directly read-off 
from the expressions below. 

Before we proceed further, we remark here that the 
above single band argument only takes into account 
bound states formed by superposition of particle and 
holes of the same band, with energies close to the Fermi 
level. For the TI in its normal phase however, there are 
bound states at E — ±vku formed by superposition of 
the conduction and valence band. Hence, our single band 
approximation for superconducting bound states is appli- 
cable only when the states at E = ±vk\\ are sufficiently 
far away from /i so that we can ignore the hydridization 
of our states with these due to the superconducting pair- 
ing A, i.e., we need — \vk\\\ 3> A. For \fi\ ^ m, we 
thus we expect that we can capture the superconducting 
bound states only when k^kp\\, the Fermi momentum. 
More discussion on this will be given below. 

Now we apply our above results to the odd-parity 
phases. 

Ai u : d(k) is available in eq and (|4"0"1) . We write 
k\n k x x ~\~ kyij ~\~ k z z and fcout = k x x -t- kyy k z z , k z ^> 0. 
Since the system is rotationally symmetric about z, let 
us consider k y — 0. We get 

/ \ t* 7* ~ i a i9 7Y17J zk^vkx ^ / , r,\ 

(sgn/x)d in x d out = 2\A lu \ z y (48 

(m 2 + vjkj)\n\ 

with vn x k — vk x y. We see that the dispersion is 
RN if sgn(m?; z ) > 0, but AN if sgn(mv z ) < 0. The 
magnitude of the group velocity of the bound state, 
\dEt/dk\\ |, can also be obtained easily. For small k x , 

sinC = d in x <i out /|A'| 2 w (tt — (), since C is close to n, so 
\E b \ a |A'|(tt - 0/2. Using |A'| = A lu (l - (m/^) 2 ) 1 / 2 , 
we get 

TOH^iM^^y" 2 * (49) 



d(k) = A u { K ri f-z+ sgn M 



The factor v 2 (k 2 — k 2 )/[j, 2 only gives a small correction 
to the results below and will be ignored for simplicity. 
For the k x -k z plane, k y = 0, 

(agnfi)d in x d out = 2 A 2 mVz ^vk x „ ^ 
whereas vnxk — vk x y. Thus we have RN if sgn(mw^) > 



For |m//i| 1, this group velocity is reduced compared 
with the BW phase by a factor \m/n\ (see also^). 

In our single band description, the bound state energy 
approaches the bulk gap when fc|| approaches the Fermi 
momentum kp\\, here given by (/i 2 — m 2 ) 1 / 2 /|t>|, as the 

effective phase difference x vanishes when k[ n and fc out 
becomes parallel. The bound state spectrum is thus of 
the form in Fig [TJa) or Fig HJb) according to whether 
sgn(muj)^0 (when v taken as < according to SecHI)). 
In contrast, the bound state spectrum in the full two- 
band calculation for sgn(mw z ) < is schematically shown 
in Fig [T] (c) 21 i 22 . Though our one-band model captures 
correctly the sign of the group velocity for feii^/cpM, it 
does not capture the behavior at larger fcii . This is in 
retrospect not surprising. For large , the Cooper pair- 
ing plays no role, and the sign for the dispersion must 
be the same as the corresponding normal phase. In the 
later case, when the system is a TI, sgn(rTi^) < and 
the positive energy branch has spin along vz x k. Hence, 
for sgn(mu 2 ) < 0, the sign of the dispersion at small 
and large fey must be opposite, as in FigQTc). This sign 
change has been found earlier by other authors 2 ^— . An 
alternative view of this sign change has been given by-2 2 .. 
(see also Appendix [Al for more discussions.) 

Am'- d{k) is given in eqn (|43[) . which is independent of 
the sign of m and v z . Since d(k) is independent of k Zl 
dm = d out . There are no surface bound states (in our 
approximation where the order parameter is constant up 
to the surface). This conclusion is in agreement with Fig 
5 a, c ofSi. Note that the gap magnitude vanishes for 
normal incidence, where k x = k y = 0. 

E u . We consider the state f, 2 t> — 11 1,2 |>), the 
first line under E u in Table U (The result for the second 
line is the same except for a tt/2 rotation about z). d(k) 
is available in cq (|46|) . For near normal incidence, it can 
be rewritten as 



- i / \ v zk z v 2 k x k y „ I 

x + {sg ^2UA^^ y j (50) 
I 

and AN if vice versa. For the later case, as argued in the 
last paragraph for A\ u , the dispersion should be the same 
as the normal phase for large fc|| , hence we again expect 
a sign change for the group velocity at some k\\ . This is 
in accordance with Fig 6a of^L (recall^!) and o^. 

Before we depart from this section, we would like to 
make a remark on the more general case where the mo- 



v z k z 



1 v 2 {k 2 x 
1 + - 



11 



mentum dependence of the term m in eq ([I]) is included. 
All the above single-band calculations can be simply gen- 
eralized to this case. For weak superconducting pair- 
ing, it is the value of m(k) at the Fermi surface that is 
physically relevant. For a TI, m(k) — mo + Ck 2 where 
mo and C has opposite signs (We shall ignore possible 
anisotropies in the C term as they do not affect the argu- 
ments below). Thus it is possible that the sign of m(k) for 
k at the Fermi surface be different from m{k = 0) = mo. 
Hence, sgn(m(fc)i> 2 ) can be positive for a TI, even though 
sgn(mof z ) < is required, and the spectra can change 
from AN to RN with increasing or \C\ (for Ai u and 

2 

E u ), at the point where mo + = in the ^> \m\ 

limit. This is a simple explanation of the finding of2£. In 
this respect, thus rigorously speaking, RN or AN is not 
an indication of topological character of the underlying 
normal phase, but rather the d(k) configuration of the 
superconducting phase. 

V. CONCLUSION 

In this paper, we have constructed a pseudospin basis 
to describe the normal state of Bi2Se3. Superconductiv- 
ity is then expressed in this basis. Using this approach, 
many of our previous knowledge in unconventional super- 
conductivity can then be directly applied, especially for 
the bulk. We have also shown that many features of the 
surface bound states can also be understood in this way. 
Although we have concentrated on the surfaces parallel 
to the Bi2Se3 quintuple layers in this paper, the same 
considerations are applicable to other surfaces as well as 
systems with other symmetries. This picture however 
misses the topological properties of the normal phase, 
which in turn some features of the surface bound states in 
the superconducting phases, but only at momenta com- 
parable with or larger than the Fermi momenta in the 
weak-coupling limit. 
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Appendix A: 

We consider discuss some topological aspects and the 
surface states for the A iu and A ig phases. Let us first 
consider the Ai u phase, and begin with simple conti- 
nuity arguments. For simplicity, we shall consider non- 
vanishing v and v Zl varying only m. For a single-band 



model with d{k) in eq (|39f and (|40|) . we note that at 
m = 0, d(k) only has z component and is odd in k z . 
The bulk state then has a line node on the equator, and 
the surface state is simply a flat band independent of k\\ . 
This is how the surface state spectra evolve between Fig 
QJa) and FigQJb) when m changes sign. 

Time-reversal symmetric superconductors in three- 
dimension can be characterized by a winding number 
44-46 xhis value can be evaluated by first transforming 
the Hamiltonian into an off-diagonal form. The Hamilto- 
nian for the A\ u phase in the Nambu-II notation is, in the 
one-band model, £t z + (d(k) ■ s)t x . Here £ is the kinetic 
energy measured with respect to the chemical potential 
(i. For example, for a quadratic band with particle-like 

dispersion, £ — — [i with /i > 0, whereas for a hole- 
,2 

like band, £ = —^m ~ f 1 w ^ tn M < 0- Here M > is an 
effective mass (not to be confused with m in eq ([T])). The 
Hamiltonian becomes off-diagonal under a rotation in r 
space, such as r z — > t x , t u 
becomes 



T y . Then the Hamiltonian 



H s = 





4 



h s 




(Al) 



where hg = £ + i(d(k) ■ s). The winding number can be 
evaluated from 



W = 



1 



24tt 2 



d 3 ke abc Tr 



.dq_ ^dq_ t dq 
dk a dkb dk c 



(A2) 



where a, b,c = x,y,z, e a i, c is the fully antisymmetric 
tensor, and q = hs/Eg is a unitary matrix. Here 
E s = [£ 2 + d ■ d] 1 / 2 . For the BW phase with d paral- 
lel to k, we get W — sgn^j£2 For our state with eq 
and 1|40|). we get instead 



W = (sgn/i)sgn(A lu w z ) 



(A3) 



independent of the sign for m, since both d x , y would 
change sign under a sign change of m. Thus, the wind- 
ing number seems insufficient to indicate the possible 
change in the surface state spectra (and the associated 
bulk topology) between Fig [1] (a) and (b) . 

Now we turn to the full two-band model, and again 
first employ only continuity arguments. The state is 
gapped if A lu ^ so long as [i 2 + A 2 U > m 2 . Hence, one 
can change sign of m without going through any gapless 
phase, provided the above inequality is satisfied. Hence 
when sgn{mv z ) changes sign, it cannot affect which spin 
species is connected to the Es > band at large kn, 
even though the sign of the dispersion can change for 
smaller fey . Hence, continuity argument shows that Fig 
[TJc) should evolve to Fig [TJa) when sgn(mv z ) changes 
from < to > 0. 

We can also examine the winding number W in the 
full two-band model. The Hamiltonian in the Nambu-II 
notation Hs = {Hn — (x)t z + Ai u a y s z T x becomes off- 
diagonal by the same rotation in f space mentioned be- 
fore, with hs = (-ffjv — m) + i^iu&yS z - Since the state 
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is gapped so long as /i 2 + A 2 U > m 2 and the winding 
number cannot change within a gapped phase, one can 
first turn off m and then fx, provided A.± u =^ 0. The 
eigenvalues in eq (|42p become degenerate and we can use 
q = hs/Es, with q unitary, in eq (|A2[) . A direct evalua- 
tion gives W = 1 (for Ai„ > 0) irrespective of the sign 
of sgn(mi>2). There is no topological phase transition 
between sgn(mt> z ) < and > within the A± u phase. 
We note that the single band model therefore produces 
a spurious topological change when m changes sign. 

Lastly, we consider the even parity A\ g phase. In the 
single-band picture, this is the ordinary s-wave super- 
conductor. There is no phase difference for the order 
parameter at k- ln and fc out , and so no bound states are ex- 
pected. The surface states for A' lg and A lg phases have 
been investigated by^i. They showed that there are no 
surface states for A' lg , in accordance with above. Inter- 
estingly, they found that surface states survives for A lg . 
The surface states seem to be in the form of two Dirac 



cones, related by the particle-hole symmetry of the su- 
perconductivity, and crossing each other at the chemical 
potential. We do not have a simple explanation of this 
in our single-band picture. Hao and Lee^i noted that 
the A x term docs not split the crossing at the Fermi 
level. It is unclear whether this reflects any topology of 
the A lg phase, such as a possible even winding number 
that is allowed for a time-reversal symmetric even parity 
superconducting stated. We here simply note that, if 
a general A\ g is considered, the pure A lg state and the 
pure A^ g state are connected in the sense that one can 
find a path in parameter space which connects them with 
the state remaining fully gapped (see Sec lIV Ap . The A' lg 
phase is topologically trivial, and no surface states are 
expected. Thus the surface states of A x are expected 
to be destroyed in general once a finite A^ is introduced. 
This is in accordance with the argument o& where they 
showed that A^ would introduce a finite matrix element 
coupling the two surface states of A 1 



1 L. Fu, C. L. Kane and E. J. Mele, Phys. Rev. Lett. 98, 
106803 (2007) 

2 L. Fu and C. L. Kane, Phys. Rev. B 76, 045302 (2007) 

3 Y. Xia, D. Qian, D. Hsieh, L. Wray, A. Pal, H. Lin, A. 
Bansil, D. Grauer, Y. S. Hor, R. J. Cava, and M. Z. Hasan, 
Nature Phys. 5, 398 (2009). 

4 H. Zhang, C.-X. Liu, X.-L. Qi, X. Dai, Z. Fang and S.-C. 
Zhang, Nature Physics 5, 438 (2009) 

5 D. Hsieh, D. Qian, L. Wray, Y. Xia, Y. S. Hor, R. J. Cava 
and M. Z. Hasan, Nature (London) 452, 970 (2008) 

6 D. Hsieh, Y. Xia, L. Wray, D. Qian, A. Pal, J. H. Dil, J. 
Osterwalder, F. Meier, G. Bihlmayer, C. L. Kane, Y. S. 
Hor, R. J. Cava and M. Z. Hasan, Science 323, 919 (2009) 

7 Y. L. Chen, J. G. Analytis, J.-H. Chu, Z. K. Liu, S.-K. Mo, 
X. L. Qi, H. J. Zhang, D. H. Lu, X. Dai, Z. Fang, S. C. 
Zhang, I. R. Fisher, Z. Hussain, Z.-X. Shen, Science, 325, 
178 (2009) 

8 D. Hsieh, Y. Xia, D. Qian, L. Wray, J. H. Dil, F. Meier, J. 
Osterwalder, L. Patthey, J. G. Checkelsky, N. P.Ong, A. 
V. Fedorov, H. Lin, A. Bansil, D. Grauer, Y. S. Hor, R. J. 
Cava and M. Z. Hasan, Nature 460, 1101 (2009) 

9 X.-L. Qi, T. L. Hughes, S.-C. Zhang, Phys. Rev. 78, 195424 
(2008) 

10 X.-L. Qi, R. Li, J. Zang and S.-C. Zhang, Science, 323, 
1184 (2009) 

11 Y. S. Hor, A. J. Williams, J. G. Checkelsky, P. Roushan, 
J. Seo, Q. Xu, H. W. Zandbergen, A. Yazdani, N. P. Ong, 
and R. J. Cava, Phys. Rev. Lett. 104, 057001 (2010) 

12 L. A. Wray, S.-Y. Xu, Y. Xia, Y. S. Hor, D. Qian, A. V. 
Fedorov, H. Lin, A. Bansil, R. J. Cava, and M. Z. Hasan, 
Nat. Phys. 6, 855 (2010). 

13 M. Kriener, K. Segawa, Z. Ren, S. Sasaki, and Y. Ando, 
Phys. Rev. Lett. 106, 127004 (2011) 

14 Bi2Te3 is also found to be superconducting under 
pressure 1 ^. However, the bulk carriers are holes which are 
located away from the k = L point, a rather different 
situation from Cu:Be2Se3. We shall not discuss this com- 
pound here any further in this paper. 



15 J. L. Zhang, S. J. Zhang, H. M. Weng, W. Zhang, L. X. 
Yang, Q. Q. Liu, S. M. Feng, X. C. Wang, R. C. Yu, L. Z. 
Cao, L. Wang, W. G. Yang, H. Z. Liu, W. Y. Zhao, S. c. 
Zhang, X. Dai, Z. Fang and C. Q. Jin, Proc. Acad. Sci. U. 
S. A. 108, 24 (2011) 

16 S. Sasaki, M. Kriener, K. Segawa, K. Yada, Y. Tanaka, M. 
Sato, and Y. Ando, Phys. Rev. Lett. 107, 217001 (2011). 

17 T. Kirzhner, E. Lahoud, K. Chaska, Z. Salman, and A. 
Kanigel, B86, 064517 (2012) 

18 T. V. Bay, T. Naka, Y. K. Huang, H. Luigjes, M. S. Golden, 
and A. de Visser Phys. Rev. Lett. 108, 057001 (2012) 

19 N. Levy, T. Zhang, J. Ha, Fred Sharifi, A. A. Talin, Y. 
Kuk and J. A. Stroscio. larXiv:1211. 02671 

20 L. Fu and E. Berg, Phys. Rev. Lett. 105, 097001 (2010) 

21 L. Hao and T. K. Lee, Phys. Rev. B 83, 134516 (2011) 

22 T. H. Hsieh and L. Fu, Phys. Rev. Lett. 108, 107005 (2012). 

23 A. Yamakage, K. Yada, M. Sato and Y. Tanaka, Phys. 
Rev. B 85 180509(R) (2012) 

24 P. W. Anderson, Phys. Rev. B 30, 1549 (1984) 

25 G. E. Volovik and L. P. Gorkov, Pis'ma Zh. Eksp. Teor. 
Fiz. 39, 550 (1984) [JETP Lett. 39, 12 (1984)]; Zh. Eksp. 
Teor. Fiz. 88, 1412 (1985) [Sov. Phys JETP 61, 843 (1985)]. 

26 K. Ueda and T. M. Rice, Phys. Rev. B 31, 7114 (1985). 

27 E. I. Blount, Phys. Rev. B 32, 2935 (1985). 

28 M. Sigrist and K. Ueda, Rev. Mod. Phys. 63, 239 (1991) 

29 S. Yip and A. Garg, Phys. Rev. B 48, 3304 (1993) 

30 Robert Joynt and Louis Taillefer, Rev. Mod. Phys. 74, 235 
(2002) 

31 Different models have also been employed in the litera- 
ture, e.g.,— i 16 ' 21 , hence some modifications must be taken 
to compare those results with here. In Model I o^-, the 
spin part in eq ^ was instead written as k x s x + k y s y . 
Hence all the spins ir>21 must be rotated by ir/2 about z in 
order to compare with ours here. Fortunately, this would 
not alter anything for pairings between opposite spins, that 
is, for the Ai g , Ai u and A2 U phases in Table[l] For the two- 
dimensional odd parity state E u which involves pairing be- 
tween parallel spins, this rotation brings in the factors =fz 



13 



for |1 t 2 t> an d 1 4- 2 i> respectively, thus interchanges 
the two entries under E u in Table U We shall not try to 
provide the comparison with other models in this paper. 
C.-X. Liu, X.-L. Qi, H. J. Zhang, X. Dai, Z. Fang and S.-C. 
Zhang, Phys. Rev. B 82, 045122 (2010) 
F. Zhang, C. L. Kane and E. J. Mele, Phys. Rev. B 86, 
081303 (2012) 

We are using this term with a different meaning from 
the usual "two-band superconductivity" in the literature. 
Here, we only have one Fermi surface for the corresponding 
normal phase in the absence of pairing potentials. 
A. J. Leggett, Rev. Mod. Phys. 47, 331 (1975) 
There are no entries listed under columns (iii), (iv) and 
(v) for A2 g and E g , since they do not arise if we insist 
that there is no momentum dependence when the pairs are 
expressed in the orbital and spin basis. They can however 
occur if we allow more general momentum dependence. 
One is free to work in other basis, e.g.—, but the physical 
picture needed would be different. Also, one would not be 
able to employ the tables in22r—. 

K. Michaeli and L. Fu, Phys. Rev. Lett. 109, 187003 (2012) 
We shall not worry about how we define states at k z — 0. 
We can define these by approprate limiting procedures. 
This can lead to discontinuities of these wavefunctions as 
functions of k (see^) , but these complications would not 
occur for the pair wavefunctions, our major concern in this 
paper. 

The pseudospin pj obeys Pj(k) — > pj(—k) under parity 
while pj(k) —pj(—k) under time-reversal. This is appli- 
cable to both pi, 2, 3 and px, y ,z defined below due to eq ((2} 
and ©. 

We are using 1, 2, 3 here for these matrices in primed basis 
since they do not have the rotational properties of x,y,z. 
R. Balian and N. R. Werthamer, Phys. Rev. 131, 1553 
(1963) 

N. P. Butch, P. Syers, K. Kirshenbaum, A. P. Hope and J. 

Paglione, Phys. Rev. B 84, 220504(R) (2011) 

A. P. Schnyder, S. Ryu, A. Furusaki, and A. W. W. Lud- 

wig, Phys. Rev. B 78, 195125 (2008). 

M. Sato, Phys. Rev. B 79, 214526 (2009) 

X.-L. Qi, T. L. Hughes, S.-C. Zhang, Phys. Rev. 81, 134508 

(2010) 

S.-K. Yip J. Low Temp. Phys. 160, 12 (2010) 
We can separate the reflection due to the vaccum-sample 
surface and the Andreev reflection due to the order param- 
eter variations by imagining inserting a thin normal layer 
near the surface, as done in, e.g.— . Since the reflection 
coefficient for the hole is given by the complex conjugate 
of that of the particle, the phase shift from the reflec- 
tion coefficient iig of Sec Mil cancels, if we ignore the small 
wavevector difference between the particle and hole. 
C. W. J. Beenakker, Phys. Rev. B 46, 12841 (1992) 
These results can be derived easily from solving the An- 
dreev equation, or locating the poles of the Green's func- 



tion in, e.g.,—. 

I. O. Kulik and A. N. Omelyanchuk, Sov. J. Low Temp. 3, 
459 (1977) 

This sign change arises from Ek ~ E^f + {dE^ / dk)(k — 
kkF), when one derives the Andreev equation. 
There seem to have an ambiguity in the sign of W since 
one can characterize the same phase with d antiparallel to 
k by making use of a gauge transformation. 




FIG. 1: Schematic dispersions for surface bound states in 
the Aiu phase, (a): Single-band or full two-band picture, 
sgn(mw z ) > 0, (b): Single-band picture, sgn(mti z ) < 0, (c): 
full two-band picture, sgn(rm; z ) < 0. Here, the arrow head 
(tail) indicates that the spin is pointing out of (in to) the 
plane, v < 0. 



